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Implications  of  conservation  equations  for  the  determination 

of  absolute  velocities 

by  C.  L.  Frenzen1  and  Alex  C.  Wara-Vamas2 


ABSTRACT 

The  consequences  of  assuming  that  density  is  conserved  in  the  problem  of  determining 
absolute  velocities  are  investigated.  Two  questions  are  considered:  (i)  the  constraints  that  the 
density  must  satisfy  to  be  compatible  with  the  assumed  geostrophic  and  hydrostatic  dynamics 
and  (ii)  whether  and  to  what  extent  the  indeterminacy  in  this  dynamics  is  removed  by  this 
additional  assumption. 


1.  Introduction 

With  the  dawn  of  ocean  forecasting  approaching,  the  need  to  know  velocity  fields  for 
use  as  initial  conditions  in  a  prediction  scheme  is  keenly  felt.  But  direct  velocity 
measurements  are  still  difficult  to  make,  and  so  we  are  forced  to  rely  on  the  classical 
procedure  which  consists  of  inferring  the  velocity  fields  from  temperature  and  salinity 
measurements.  In  this  procedure,  we  are  given  the  density  field  p(x,  y,  z),  and  we  must 
find  velocity  fields  u(x,  y,  z),  v(x\  y,  z ),  w(x,  y,  z)  and  a  pressure  field  p(x,  y,  z)  such 
that 


fv  =  Po'p,  (11a) 

fu^-Po'p,  (1.1b) 

0  -  -P~'Pz  ~  g  (1.1c) 

ux  +  +  =  0  ( 1 .  Id) 


The  notations  in  the  above  equations  are  standard:  x,  y,  z  are  cartesian  coordinates;/ 
is  the  Coriolis  parameter  and  g  is  the  gravitational  acceleration;  p0  is  an  average  value 
of  the  density. 

As  is  very  well  known,  the  integration  of  the  hydrostatic  equation  from  the  ocean’s 
bottom  at  z  -  -  h  to  a  depth  z ,  viz. 

P  -  -g  ['  p(x,y,z’)dz'  +  v(x,y)  (1.2) 
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introduces  an  unknown  function  of  x,  y  denoted  here  by  ir.  This  ‘barotropic’  component 
of  the  pressure  cannot  be  determined  from  the  dynamics  (1  .la,  b,  c,  d). 

The  indeterminacy  in  the  pressure  entails  a  similar  indeterminacy  in  the  velocity 
fields.  Indeed,  if  we  define 


then 


where 


and 


P(x,y,z)  =  -g  f1  p{x,y,z')dz'  (1.3) 


u(x,y,z)  =  u'(x,y,  z)  +  u(x,y) 

v(*,  y,  z)  =  v'(x,  y,  z)  +  v(x,  y) 

(1.4) 

u'  =  -(Pof)-'Py 

(1.5a) 

v'  =  (pjy'px 

(1.5b) 

u - (Pofr'lTy 

(1.6a) 

v  -  (po/y'*x 

(1.6b) 

The  vertical  velocity,  which  is  obtained  by  integrating  the  continuity  equation  in  the 
vertical. 


w  =  w(x,y, -h)  -  (ux  +  vy)  dz’  (1.7) 

J -h 

also  contains  the  same  indeterminacy.  The  search  for  the  famous  ‘level  of  no  motion’  is 
intimately  related  to  the  determination  of  v{x,  y). 

Following  in  the  footsteps  of  Worthington  (1976),  Wunsch  (1978),  Fiadeiro  and 
Veronis  (1982)  and  others,  we  want  to  examine  whether  the  addition  of  extra 
information  in  the  form  of  the  conservation  of  the  density  field,  viz. 

upx  +  vpy  +  wPz  =  0  (1.8) 

can  determine  x(jc,  y). 

Mathematically,  it  is  not  clear  whether  Eq.  (1.8)  is  compatible  with  the  dynamics 
(1.1).  We  shall  examine  this  question  first  on  the  /-plane  and  then  on  the  /8-plane.  We 
shall  show  that  in  both  cases  the  density  field  must  satisfy  certain  constraints.  If  these 
constraints  are  satisfied,  then  for  the  /3-plane  case,  but  not  for  the  /-plane,  t  can  be 
determined.  Our  analysis  is  carried  out  in  the  simplest  case  of  a  flat  bottom  ocean. 
However,  it  can  be  extended  to  the  case  where  bottom  topography  exists  as  well  as  to 
cases  where  quantities  other  than  density  are  conserved. 
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2.  A  mathematical  result 

An  essential  step  in  our  analysis  consists  in  finding  a  z-independent  solution  of  an 
equation  of  the  form 

A  ( x ,  y ,  z)vx  +  B (x,  y,  z)vy  =  1  (2.1) 

which  holds  in  a  cylindrical  basin  D  =  {(x,  >>,  z) :  (x,  y)  £:  D,  -h  <  z  <  0},  where  D  is 
the  horizontal  cross-section.  If  A  and  B  themselves  are  independent  of  z,  then  the 
existence  of  such  solutions  is  not  too  surprising.  Of  greater  interest  is  the  case  in  which 
such  solutions  exist  even  though  A  and  B  are  z-dependent,  i.e.  when 

Al  +  Bl^OinD.  (2.2) 

Relabeling  coordinates  if  need  be,  we  interpret  (2.2)  to  mean  that 

= M  in  D.  (2.2') 

With  the  above  conditions  we  have  the  following  result. 


a.  Theorem.  There  exists  a  z-independent  solution  w(x,  y)  of  (2.1)  if  and  only  if  for  all 
(*.  y,z)ED 


(i)  A  (x,  y ,  z)  -  k0(x ,  y)B(x,  y,  z)  +  /0(x,  y) 

(2.3a) 

(»0  /o(x,y)  *  0 

(2.3b) 

(iii)  dyUo  ')  ~  dAko/lo)  -  0. 

(2.3c) 

Furthermore,  this  solution  x(x,  y)  satisfies 

rx  =  lo'(x,y) 

(2.4a) 

*y  =  -k0(x,y)/lAx,y). 

(2.4b) 

b.  Proof.  Suppose  (2.3a,  b,  c)  hold  in  D.  We  must  show  that  a  z-independent  solution 
ir(x,  y)  of  (2.1)  exists,  and  that  it  satisfies  (2.4a,  b). 

In  view  of  (2.3b),  we  can  write  (2.3a)  thus: 


Alt'  -  Bkf0'  =  1.  (2.5) 

Now,  (2.3c)  implies  the  existence  of  a  function,  say  r(x,y),  such  that  (2.4a,  b)  hold. 
Substituting  / o'  and  -k0l o'  by  their  expressions  in  terms  of  r  in  (2.4),  we  conclude 
that  this  function  r  is  a  solution  of  (2.1). 

Conversely,  suppose  that  it  is  a  z-independent  solution  of  (2.1)  and  that  (2.2')  holds. 
Then  we  can  differentiate  (2.1)  with  respect  to  z  and  write 


Ax*x  +  B,ry  =  0. 


(2.6) 
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We  note  at  this  stage  that  AZB  -  B,A  ^  0  in  D.  Indeed,  if  this  were  the  case,  then 

B2(Aitx  +  Bvy)  -  B2, 

which  is  obtained  by  multiplying  (2.1)  by  B2,  could  be  written  as 

B(Azirx  +  B2iry)  =  B2. 

But  as  a  result  of  (2.6),  we  reach  a  contradiction.  Therefore 


a2b  -  b2a±  0 

(2.7) 

and  as  a  result,  we  can  solve  (2. 1 )  and  (2.6)  for  irx  and  iry  to  find 

b2 

ab2  -  ba2 

(2.8a) 

A. 

Xy  ab2  -  ba2 

(2.8b) 

This  step  is  not  unlike  that  taken  by  Stommel  and  Schott  ( 1 977)  and  Needier  ( 1 985)  in 
their  study  of  this  problem. 

We  must  next  insure  that  these  expressions  for  v  are  indeed  2-independent.  By 
considering  their  ratio,  we  immediately  see  that 

A,  -  k0(x,y)B2 


and  after  integration  we  arrive  at  the  condition  (2.3a)  between  A  and  B.  Substituting 
this  condition  in  (2.8a,  b),  we  arrive  at  (2.4a,  b).  Condition  (2.3c)  arises  from  forming 
nxy  in  two  different  ways.  Finally,  the  condition  (2.7)  translates  into  (2.3b).  We  can 
easily  verify  that  further  z  differentiations  do  not  introduce  other  conditions.  This 
completes  the  proof. 

In  closing  this  section,  we  give  an  example  of  a  partial  differential  equation  of  the 
form  (2.1)  in  which  the  coefficients  A  and  B  satisfy  all  the  conditions  of  the  theorem. 
The  example  is: 


B(x,y,  z)  + 


1 

fix  -  y) 


xx  +  B(x,y,  z)iry  =  1 


where  B  is  an  analytic  function  of  z.  In  this  case, 


k0  =  1 

, _ i _ 

0  f'(x-y) 


and 


*  -  fix  -  y). 


1988] 


Frenzen  &  Warn-  Varnas:  Conservation  equations 


705 


3.  The  /-plane 

We  consider  in  this  section  the  case  of  the  /-plane,  simply  as  a  foil  for  the  /3-plane 
which  is  discussed  next.  In  this  case,  the  flow  is  horizontally  nondivergent.  This, 
together  with  the  fact  that  the  bottom  is  flat,  implies  that 

w  =  0.  (3.1) 

Thus,  the  added  conservation  equation  reads 

«Px  +  vpy  =  0.  (3.2) 

Making  use  of  the  hydrostatic  and  geostrophic  balances,  this  equation  can  be  written  as 
follows: 


-PyPxz  +  PxPyx  ”  0  (3.3) 

and  since  p\  +  p\  is  not  identically  zero  in  D  lest  the  ocean  would  be  at  rest,  this  implies 
that 


px=  -py  tan  a0(x,y)  (3.4) 

or  equivalently 

v(x,  y,  z)  =  u(x ,  y,  z)  tan  a0(x,  y).  (3.5) 

Thus,  at  a  given  latitude  and  longitude,  the  velocity,  which  is  purely  horizontal,  has  the 
same  direction  for  all  depths!  This  result,  in  turn,  places  some  constraint  on  the  data. 
Indeed,  substituting  (3.5)  in  (3.2)  shows  that 

Px  4-  py  tan  a0  =  0.  (3.6) 

This  constraint  is  typical  of  the  constraints  the  conserved  quantity  must  satisfy  for 
the  data  to  be  compatible  with  the  assumed  dynamics.  What  this  means  in  practice  is 
that  at  each  station,  z-independent  a0  are  generated  by  using  the  best  fitting  line 
passing  through  points  with  abscissas  py  and  px  (or  whatever  the  corresponding 
quantities  are  for  the  case  at  hand). 

As  we  shall  see,  a  similar  constraint  arises  on  the  /3-plane.  Having  satisfied  this 
constraint,  are  we  now  able  to  determine  ir?  The  answer  is  unfortunately  no.  Indeed, 
from  the  definition  ( 1 .3)  of  P,  it  follows  that: 

Px  +  Py  tan  a0  =  0.  (3.7) 

This  means  that  the  baroclinic  components  of  the  flow  field  satisfy  the  conservation 


706  Journal  of  Marine  Research  [46, 4 

equation  by  themselves.  Therefore,  the  need  for  a  barotropic  correction  is  not 
imperative.  From  (3.4)  and  (3.7)  it  follows  that 

xx  +  it y  tan  a0  =  0.  (3.8) 

If  we  were  forced  to  proceed  further,  then  we  would  need  to  know  the  distribution  of 
sources  and  sinks  at  the  edge  of  our  basin  in  order  to  solve  this  equation  for  x. 


4.  The  0-plane 

Because  of  the  /3-effect,  the  horizontal  mass  divergence  is  -0v//o  and  therefore  the 
vertical  velocity  as  given  by  (1.7)  is 


w(jf,  y,  z)  =  0/o  1  fz  v(x,y,  z')  dz' 

J-h 

or  in  terms  of  P  and  w 

w  =  (0/foPo)  {yj  p,  dz'  +  [z  +  h]rx 
We  shall  find  it  convenient  at  times  to  write 


(4.1) 


(4.2) 


fzPxdz'. 

-h 


(4.3) 


The  conservation  equation  for  p  now  becomes 

ip,  +  (0[*  +  h]/fo)pz)irx  -  pxry  =  -pyPx  +  pxPy  -  lp0fo)pxW.  (4.4) 

The  question  once  again  is  whether  the  right-hand  side  of  (4.4)  vanishes,  i.e.  whether 
the  baroclinic  fields  identically  satisfy  the  conservatiort  equation.  We  recall  that  for 
d  =  0  the  answer  to  this  question  was  always  yes.  If  each  side  of  (4.4)  is  separately  zero, 
then  the  density  field  which  makes  up  the  data  is  constrained  to  satisfy  two  distinct 
equations.  The  equation  stemming  from  the  baroclinic  part  is 

-PyPx  +  pxPy  -  (Po/o)Pr^  -  0  (4.5) 

whereas  that  from  the  barotropic  part  yields 

Px  +  i Py  +  (01*  +  A]//oW  tan  oo  =  0  (4.6) 


where  a0  is,  as  previously,  solely  a  function  of  x  and  y.  Are  there  isopycnal  surfaces 
which  satisfy  both  of  the  above  equations? 

We  shall  show  that  this  question  must  be  answered  negatively.  Integrating  (4.6)  over 
z  and  using  the  definition  (1.3)  of  P,  we  deduce  that 


/\  + 


0(z  +  h) 

fo 


Pg+jP 

Jo 


0. 


Px  +  tan  a 0 


(4.7) 
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Next,  we  eliminate  py  and  Py  from  (4.5)  using  (4.6)  and  (4.7).  The  result  is: 


Pfo'px 


p  +  pg  jf  *  dz'  =  0/0  y  -Px /'  dz'  +  jf  *  Px  dz' 


If  we  differentiate  this  equation  with  respect  to  z,  we  see  that 


fifo  'Pxx  l 

\p  +  pg  f*  dz'\ 

=  Bfo'Pzx  ' 

\-Px  f^z’  +  f’Pxdz’ 

[  J-k  j 

J-„ 

The  above  two  equations  imply  that 


i.e. 


PzPxz  =  PxPzz 


Px  =  px  tan70(x,>-). 
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(4.8) 


(4.9) 


(4.10) 


We  have  thus  reduced  the  problem  to  finding  whether  there  exists  a  function  p(x,  y,  z) 
which  satisfies  the  system  of  equations 


£,  [p]  *  tan  70 px  -  p2  -  0 

£ i  [p]  =  Px  +  tan  a0py  +  (0[z  +  h]/f0)  tan  a0pz  -  0. 

At  this  stage  we  appeal  to  the  theory  of  systems  of  first  order  partial  differential 
equations  (see  e.g.  Smirnov,  1964,  p.  34.).  The  first  step  in  the  theory  consists  in 
attempting  to  derive  further  independent  equations  for  p.  Such  additional  equations 
are  obtained  by  forming  Poisson  brackets.  Thus,  denoting  temporarily  x ,  y,  z  by  xx ,  x2, 
x3  and  px,  py,  pz  by  pt,  p2,  p},  the  theory  leads  us  to  consider 

£3[p]  -  E  {(dE,/dXl)(dE2/dPi)  -  0£2/dx,)(d£,/d/>,)}  -  0 

I- I 


i.e. 

-  [(d-Yo/dx)  +  (dy0/dy)  tan  a0(cos  7o)"2]Px  +  [( da0/dy )  tan  y0(cos  ao)-2]p,, 

+  [f)(z  +  h)f~\da0/dy)  tan  70(cos  oq)-2  -  /3/_l  tan  «o]pr  =  0  (4.1 1) 

This  equation  is  independent  of  the  other  two.  Having  obtained  three  independent, 
homogeneous,  first  order  partial  differential  equation  for  p  which  is  a  function  of  three 
independent  variables,  the  theory  tells  us  that  p  =  constant  is  the  only  possible 
solution. 

Thus,  on  the  (8-plane  the  baroclinic  velocity  fields  cannot  satisfy  the  density 
conservation.  Therefore,  we  can  write  (4.4)  as 


Airx  +  Bvy  =  1 


(4.12) 
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where 


and 


P,  +  P,[P(Z  +  h)/f0) 
PxPy  -  PyPz  ~  {Pof0)P;W' 


(4.13a) 


B  =  - 


PxPy 


Px _ 

PyPx  -  (Pofo)Pz  W 


(4.13b) 


If  we  assume  that  A  and  B  are  functions  of  depth,  then  the  theorem  yields  the  following 
constraint  on  data 


A  =  k0(x,y)B  +  l0(x,y ) 

and  the  usual  expressions  for  the  barotropic  components,  namely 

*x  -  /o  '(  v,  y) 

ir,  -  ~k0(x,y)/l0(x,y). 

On  the  0-plane,  the  density  conservation  is  stringent  enough  to  determine  the 
barotropic  field  completely. 

Finally,  we  note  that  the  same  constraints  which  beta  spiral  calculations  place  on 
density  and  potential  vorticity  are  present  in  our  analysis:  potential  vorticity  cannot  be 
a  function  of  density  alone,  i.e.  material  surfaces  of  density  and  potential  vorticity 
cannot  coincide.  From  (1.1)  and  (1.8)  it  follows  that  potential  vorticity  is  conserved 
(Ertel’s  Theorem): 


uqx  +  vqy  +  wq,  =  0  (4.14) 

where  q  -  fpt.  If  (4.14)  is  used  to  eliminate  iv  from  the  density  conservation  equation 
one  obtains 


*x(P:yPz  +  ( 0/f)P x  -  P.-zPy)  -  Vy(PzxPl  ~  PzzPx)  = 

~  PxiPzyPz  +  ( &/f)pl  -  PzzPy)  +  Py(PzxP ~  PzzPx)-  (4.15) 

When  q  =  G(p)  for  some  function  G  both  sides  of  (4.15)  are  zero  and  the  barotropic 
field  cannot  be  determined.  The  same  degeneracy  is  also  present  in  (4.4):  if  q  =  G(p) 
the  right-hand  side  of  (4.4)  is  zero  and  the  barotropic  field  cannot  be  completely 
determined  from  density  conservation  alone.  One  way  to  see  this  is  to  note  that  the 
relationship 


implies  that 


fPz  =  G{p) 


(4.16) 


p(x,y,z)  =  K(z/f  +  Q(x,y)) 


(4.17) 
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Figure  1 .  Normalized  error,  r2  -(E/S)2,  in  REX  area.  E  is  the  error  of  the  estimated  field  from 
true  field.  S  is  the  signal  variance. 


where  K'(x)  -  G(K(x))  and  Q  is  some  arbitrary  function  determined  by  boundary 
conditions.  Substituting  (4.17)  into  the  right-hand  side  of  (4.4)  and  integrating  the 
resulting  expression  for  W  by  parts  then  results  in  cancellation  of  all  terms  present 
there.  Consequently  the  barotropic  field  cannot  be  completely  determined  by  density 
conservation  alone  when  material  surfaces  of  density  and  potential  vorticity  coincide. 
We  note  that  perhaps  the  discrepancies  in  the  examples  which  follow,  one  including 
upper  ocean  data  and  the  other  not,  suggest  that  this  “ q  =  G(p)"  degeneracy  occurs  in 
the  upper  ocean  at  the  site  of  the  example.  For  further  information  about  the 
occurrence  of  this  possibility,  we  mention  the  work  of  Rhines  and  Young  (1982). 
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Figure  2.  N-S  density  gradient  vs.  E-W  density  gradient  at  37N  and  64E  for  different  depths. 
The  points  in  the  first  quadrant  represent  measurements  taken  on  August  13-15.  1 985  as  part 
of  the  REX  experiment  in  the  upper  700  m.  Points  in  the  fourth  quadrant  are  from 
climatological  data  for  depths  greater  than  700  m. 

5.  Concluding  remarks 

We  conclude  by  illustrating  the  procedure  with  the  analysis  of  a  particular  data  set. 
The  set  selected  was  obtained  in  the  REX  (Regional  Energetic  Experiment)  experi¬ 
ment.  AXBT  profiles  were  taken  on  August  13,  14,  and  15,  1985.  The  drops  were  at 
20  km  intervals  along  altimeter  tracks  spaced  about  90  km  apart. 

The  AXBT  profiles  were  smoothed,  in  order  to  remove  the  finescale  variability,  and 
combined  with  GDEM  (Generalized  Digital  Environmental  Model)  below  700  m 
climatology  depths. 

The  data  was  objectively  analyzed  for  the  three-day  period  using  optimal  estimation 
theory  (Gandin,  1963).  This  technique  interpolates  mixed  observations  taken  at 
random  points  in  space  and  time  to  a  uniform  grid.  It  uses  the  statistical  properties  of 
the  field  and  minimizes  the  difference  between  analysis  and  observations. 

The  e-folding  correlation  scales  were  assumed  to  be  half  a  degree  and  the  noise  to 
signal  ratio  was  .01.  The  analyzed  temperature  field  was  three-dimensional  and 
contained  all  the  characteristic  features  of  the  area.  The  details  of  the  optimal 
estimation  analysis  for  the  REX  area  are  described  in  Bennett  and  May  (1988). 
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Figure  3.  A  and  B  fields  as  computed  from  (4.13a,  b).  The  line  passing  through  these  points 

determines  the  parameters  /0  and  k0  at  that  location. 

The  normalized  error  obtained  from  the  optimal  estimation  analysis  is  shown  in 
Figure  1.  It  ranges  from  .4  at  AXBT  location,  indicated  by  small  triangles  to  1.0  in 
regions  of  no  data.  Signal  variance  varied  from  ,4°C  at  surface  to  .05°C  below  700  m. 

We  chose  a  region  centered  at  37N  and  64E  for  our  calculations.  This  was  an  area 
that  contained  no  surface  signature  of  the  Guff  Stream  fronts  and  eddies.  Below  the 
surface  the  effects  of  the  Gulf  Stream  were  still  present  because  of  the  tilt  versus  depth 
of  the  front  (tilts  away  from  shore).  We  used  the  optimally  estimated  fields  for  our 
calculations  with  a  horizontal  differentiation  interval  of  one  degree  in  latitude  and 
longitude. 

Figure  2  shows  a  scatter  diagram  of  px  versus  py.  The  errors  in  px  and  py  range  from 
.23  x  10-5  mks  at  surface  to  .05  x  10'5  mks  at  700  m  and  below.  The  points  laying  in 
the  first  quadrant  are  from  the  REX  data  and  represent  values  from  the  upper 
700  meters;  whereas  those  in  the  fourth  quadrant  are  climatological  data  for  the  layers 
below  700  meters.  Figure  3  shows  the  corresponding  variables  A  and  B  as  evaluated 
from  (4. 1 3a,  b)  for  all  the  data  points.  The  points  located  in  the  region  of  A  -  200  mks 
and  B  =  200  mks  correspond  to  the  upper  700  m  region. 

Figure  3  shows  the  corresponding  values  of  the  variables  A  and  B  as  evaluated  from 
(4.13a,  b)  for  all  the  data  points.  The  least  square  fit  through  these  points  yielded: 

4  =  -312.19  mks 
4/k0  =  49.62  mks 

which  implies  barotropic  velocity  components  of 

u  =  43  cm/s 
v  =  -6.8  cm/s. 
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Figure  4.  The  baroclinic  components  of  the  velocity.  To  these  must  be  added  barotropic 
components  of  43  cm/s  East,  6.8  cm/s  South. 

The  “baroclinic”  velocity  is  shown  in  Figure  4,  (obtained  by  integrating  Eq.  (1.3)  from 

-h  to  z  and  using  (1.5e)  and  (1.5b)).  Incidentally,  if  the  climatological  data  alone 
were  used,  then 


/0  =  -3219.5  mks 
W*o  =  152.4  mks 

and  the  corresponding  barotropic  velocity  components  would  be 


«  =  14  cm/s 
v  =  — .7  cm/s. 
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Measurements  during  REX,  with  an  inverted  echo  sounder  and  a  “pressure”  gauge 
yields  ratios  of  barotropic  to  baroclinic  velocities  of  the  order  of  30%  (Hallock,  pers. 
comm.).  According  to  these  measurements  our  calculated  barotropic  velocities  (u  = 
43  cm/sec,  v - 6.8  cm/sec)  are  in  the  expected  range. 

All  of  our  model  equations  are  assumed  to  hold  during  our  calculations.  The  upper 
ocean,  however,  is  characterized  by  a  mixed  layer  region.  This  region  is  approximately 
50  m  deep  during  the  REX  experiment  and  is  controlled  by  atmospheric  forcing.  Its 
representation  would  require  the  introduction  of  a  mixed  layer  model  that  would 
represent  the  time  variation  of  density  and  vertical  diffusion  of  momentum  and  heat 
flux.  Since  this  region  is  only  50  m  deep,  we  do  expect  our  model  equation  to  hold  over 
most  of  the  vertical  column.  Another  deviation  from  the  model  equations  can  occur  due 
to  a  degeneracy  in  the  density  conservation  equation  when  the  potential  vorticity  is  a 
function  of  density  alone.  In  this  case  the  barotropic  field  cannot  be  completely 
determined  by  density  conservation  since  material  surfaces  of  density  and  vorticity 
coincide.  The  numerical  analysis  suggests  such  a  possibility  due  to  differences  in 
behavior  between  the  upper  and  lower  regions  of  the  ocean. 

In  summary,  we  have  attempted  to  elucidate  the  questions  of  how  a  depth-dependent 
conserved  field  can  remove  the  depth  independent  indeterminacy  of  the  velocity  fields. 
We  have  seen  that  this  attempt  leads  naturally  to  a  characterization  of  the  class  of 
acceptable  density  fields.  In  turn,  this  suggests  a  very  natural  variational  procedure, 
namely  that  the  “distance”  between  the  data  and  this  set  of  acceptable  densities  be 
minimized.  For  the  /3-plane,  there  is  a  unique  velocity  field  associated  with  this 
projection. 
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